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ABSTRACT 

In the seconds after collapse of a massive star, the newborn proto-neutron star (PNS) radiates neutrinos of 
all flavors. The absorption of electron-type neutrinos below the radius of the stalled Shockwave may drive 
explosions (the "neutrino mechanism"). Because the heating rate is proportional to the neutrino energy, flavor 
conversion of higher-energy /i and t neutrinos to electron-type neutrinos via collective neutrino oscillations 
(CiyQ) can increase the heating rate, and potentially drive explosions. We solve the steady-state boundary 
value problem of spherically-symmetric accretion between the PNS surface (r^) and the shock (rs), for the first 
time including a scheme for flavor conversion via Ci^O. For a given r^, PNS mass (M), and accretion rate (M), 
we calculate the critical neutrino luminosity above which accretion is impossible and explosion results. We 
show that CiyO decreases the critical luminosity by a factor of ~ L5, but only if the flavor conversion is fully 
completed inside rs- The effect is smaller for partial conversion. The shock radius and the physical scale for 
flavor conversion depend differently on the parameters of the problem. We quantify these dependencies and 
find that Ci^O lowers the critical luminosity substantially for small M and M, and large r^. Thus, Ci'O can be 
important for Shockwave revival if PNSs contract slowly, which may favor a stiff nuclear equation of state, and 
if progenitors reach low M at early times after collapse, which favors the lowest-mass massive stars. 
Subject headings: supernovae: general — neutrinos 



1. INTRODUCTION 

Core-collapse supernovae announce the deaths of mas- 
sive stars. The explosion develops deep in the opti- 
cafly thick stellar core and the dynamics of the neu- 
tron star in formation reveals itself only through the 
emission of neutrinos, which are scarcely detected (e.g. 
lyuksel & Beacom 2007). Computer simulations are thus 
the primary means to probe the dynamics of the super- 
nova explosions. Unfortunately, only the lowest-mass 
proge n itors explode dRa mpp & Janka 2000; Bruenn et al. 
2001'; 'Liebendorfer et al. 2001; Mezzacappa et al. 200U 
Thompson et al. 2003; Kitaura et al. 2006; Janka et al. 20081. 
In particular, it proves difficult to revive the outward move- 
ment of the Shockwave, which forms when the collapsing 
core reaches nuclear density, and which stops its progress 
due to neutrino em ission losses and other effects (e.g. 
iBurrows & Lattinierl 1985; Bruenn 1989a b). Accretion of 
matter through the stalled shock ensues and lasts for many 
dynamical times before explosion or eventual black hole for- 
mation. 

The hot proto-neutron star (PNS) cools by emission of 
neutrinos of all flavors. Because i/^, v^, v^, and Vr (here- 
after collectively denoted as do not interact with the PNS 
matter through the charged-current interactions, they decou- 
ple from matter at smaller PNS radii and higher tempera- 
tures, and thus have higher average energy than and v^. 
A fraction of the and are absorbed below the accretion 
shock, and the associated deposition of energy plays a sig- 
nificant role in the dynamics of the supernova and perhaps 
in the revival of the Shockwave (e.g. [Colgate & Wh"itg|1966l; 
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iBethe & Wilsonlll985h . Specifically, the "neutrino mecha- 
nism", as formulated bv IBurrows & Goshvl ([T993), states that 
the steady-state accretion through the shock turns into an ex- 
plosion when the neutrino luminosity of the PNS (L,y^,core) ex- 
ceeds a critical value, L""core- In lPeicha & ThompsonI (1201 Ih 
(hereafter ' Paper J) we showed that L™'coie is equivalent to 
reaching max(c|/Vesj) ~ 0.19 in the accretion flow, where 
is the sound speed and Vese is the local escape velocity. This 
"antesonic" condition is a manifestation of the inability of the 
flow to satisfy both the shock jump conditions and the Eu- 
ler equations for the accretion flow simultaneously. We also 
determined the dependence of i™'core '^he key parameters 
of the problem, including the energies of the neutrinos over 
a wide range of parameter values. Specifically, and most im- 
portantly for this paper, we found that i™'core is proportional 
to the inverse square of the and energies. 

Simulations of supernovae generally fail because the neu- 



trino luminosities in the models never reach L™'j,g^g. For suc- 
cessful explosions, either (i) 



needs to be decreased or 



As an example of the former, multi- 
dimensional effects like convection and S ASI decrease L™' 
( Yamasaki & Yamada 2005, 2006; Murphy & Burrowslifool 
.Nordha us et al. 2010.) b v making the heating more efficient 
(e.g. Herant et a lJll994j ;[BuiTows et al. 1995; Janka & Mulle] 
ll996nFrver& W arren' '2004^; jBuras et aL ^006a) or cool- 
ing less efficient (Paper J). As an example of the latter, 
L,j,_c.„re. can be en hanced by convect ion inside the PN S (e.g. 



Wilson & Mavlg ,1988- .Bruenn & Dineva 1996; Keii"etal] 
19961) . 



Most of the heating below the shock occurs due to absorp- 
tion of and on neutrons and protons, while the more en- 
ergetic escape without much interaction. Thus, because the 
luminosities in each flavor are similar (Janka 2001), 2/3 of 
the total neutrino luminosity generated by the PNS is essen- 
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tially unused. However, due to the high density of neutrinos 
in this region, self-interaction between neutrinos becomes im- 
portant and can lead to a range of phenomena called " collec- 
tive neutri no oscillations" (e.g. Pantaleone 1992; Dua n et aT] 
|2006, 2010). In particular, there is a possibility of an insta- 
bility (Dasgupta et al. 2009) that exchanges part of the high- 
energy spectra with the and spectra, which may then 
produce significantly more heating than calculations neg lect- 
ing neutrino oscillati ons. IChakrabortv et al.l (1201 lallbl) and 
iDasgupta et al.l(l201 I h have investigated the relevance of CvO 
for the shock revival in the core-collapse simulations of sev- 
eral progenitor models. 

Given that have significantly higher energy than and 
Pe, and that i™'(,ore proportional to inverse square of elec- 
tron neutrino energy (Paper I), the flavor conversion of high 
energy to and may have significant ramifications for 
supernova explosions. In this Letter, we quantify the changes 
to i™'core when the neutrinos are subject to flavor conversion 
due to collective neutrino oscillations (Ci^O). We model the 
neutrino conversion using the simplified analytical treatment 
summariz ed in Dasgup ta et al. (20 11) coupled to the code de- 
scribed in lPaper J in Section l2l we describe our calculations. 
In Section|3] we quantify the changes to L™'core ™d the shock 
radii, and compare the magnitude of the effect of Ci^O to other 
known pieces of physics. In light of our results, we assess on 
what timescale, and in what mass range of progenitors, CvO 
are likely to be important. In Section|4] we discuss and review 
our results. 

2. METHOD 



defined as 



We extend the code developed in iPaper II to include the ef- 
fect of Ci^O. We calculate the structure of the steady-state ac- 
cretion flow between the neutrinosphere at radius ri, and the 
standoff accretion shock at rs assuming spherical symmetry. 
The key parameters of the problem are the mass accretion rate 
through the shock M, the PNS mass M, its radius rj^, the 
luminosity of the PNS core L^,^ core, and the rms energies £i, of 
the neutrinos. We assume that all six neutrino flavors have the 
same luminosity, L^^ ^ok- We solve the two-point boundary 
value problem composed of the Euler equations for the den- 
sity p, velocity v, and temperature T along with the equation 
for the electron fraction Y^. The shock radius is determined 
self-consistently as an eigenvalue of the solution. We include 
heating and cooling due to the two most prominent charged- 
current interactions. Other interactions provide only a small 
change in L™'(,ore tiPaper J) . Unlike iPaper J, we assume that 
ALy/dr = for all neutrino flavors. The outer boundary con- 
ditions of the calculation are conservation of mass, momen- 
tum and energy through the shock. We assume pressureless 
free fall of iron with 1 /4 free-fall acceleration upstream of 
the shock. We also assume that ri, is the neutrinosphere of 
i/e. We compare the effects of CvO to our fiducial calcula- 
tion with electron neutrino rms energies e^, o = 13MeV and 
ep^,o = 15.5 MeV (Tho mpson et al. 2003). 

According to Han nestad et al.l (|2006|), the flavor conversion 
and effective increase of electron neutrino energies occurs 
above a synchronization radius rsync, which is defined as 



/i('-sync) = 41} 1 + 
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end 
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Here, f2 depends on the neutrino oscillation frequency and 
neutrino energy spectra and as in Das gu pta et al. (201 1), we 
choose a typical value = 50 km"^ . The quantity T-j is the 
ratio of the net lepton asymmetry in the system to the neutrino 
flux available for oscillations, and is defined as 



(3) 



where (\) is the neutrino number flux emanating from a unit 
area of the neutrinosphere. For v^. 



(4) 



and simila rly for (^^7. and 0^. . The collective potential ijl is 
defined as (lEsteban-Pretel et al.ll2007t IDasgupta et"ani20 1 11) 



fi(r) = V2GF(t)uM (y) 



2r2- 



(5) 



where Gf is the Fermi coupling constant and (t>,^.p = 4>i,^ + (pp^ + 
44>i^^ . Equations ([T]) and (|2|i are solved for rsync and rend given 
and the neutrino number fluxes. We note here that the de- 
pendence of equation (|4]l on r^ introduces an absolute scale 
to equations ([T]) and (|2|i, and therefore rsync and rend do not 
scale linearly with r^. Instead, in the limit of rsync /'"i^ ^ 1, 
the scaling is 



' sync 



oc L^/"* r 



1/2 



(6) 



for fixed neutrino energies. The same scaling holds for rend- 

The collective neutrino flavor conversion effectively in- 
creases Eij^ and to e^^ between rsync and rend- Motivated 
by more complete studies of the physical extent of the flavor 
conversion, and as a numerical expedient, we adopt the fol- 
lowing functional form for ep^(r) 



2r (r sync + r end) 
(rir end ^ sync) 



The flavor conversion is more or less complete at radius rend 



(r) = e^^fi + (e,y^ -£,.„,()) | ^ + ^ tanh 

^ ' ^^^^ 

where we use = 20 MeV ( [Thompson et aT] l2003h . We 

choose a ~ 0.679 to have the term in curly brackets equal 
to 5% at rsync and 95% at rend- We adopt an analogous pre- 
scription for Ep^ (r). 

The neutrino energy self-consistently enters not only in the 
heating, but also in the reaction rates for the calculation of 
Fe and in the inner boundary condition on optical depth (see 
IPaper II for details). 

3. RESULTS 

In the left panel of Figure [l] we present L""coie includ- 
ing Cj^O as a function of M for M = 1.2Mq and r^ = 60km 
(red dashed line) along with the fiducial calculation (black 
solid line). The critical cu rve in th e fiducial calculation is ap- 
proximately a power law ( IPaper IT) . We see that CvO lowers 
^"''coie to ^0.65 times the fiducial value forM < 0.01 M© s"' . 
For higher M, the critical curve turns upward and for M > 
0.32 Mq s~', it essentially coincides with the fiducial calcula- 
tion meaning that CvO have little effect. 
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Figure 1. Left: critical core neutrino luminosity i""coic ^ function of M for M — X.IMq and — 60km. The solid black line shows the fiducial calculation 
(no Ci'O), while the red dashed line is with Ci/O included. The vertical red dotted lines mark different regimes of the effect (see text and the right panel). The 
blue dash-dotted fine includes CvO and has the neutrino cooling reduced by a factor of 2 to approximate multi-dimensional effects. Right: effect of CuO on 
shock radii. Dots mark f j"' for calculations from Paper I, illustrating that f j"' depends primarily on i."."coic- Light blue, orange, and green dots are for ri, = 20, 
40, and 60 km, respectively. The black solid line shows " for the fiducial calculation from this paper for the same M and as in the left panel. The grey lines 
show rg as a function of Liy^^coTc for different M with CuO included. The grey lines terminate at ", which are connected with a dashed red line (compare with 
red dashed line, left panel). The black dotted lines show fsync and rend. For " < rsync, CuO has no effect; for r'^"^ > r^^i, the effect of the CuO is maximized. 



The behavior seen in the left panel of Figure[T]is non-trivial 
even in our simple setup, because rjync and rend are a function 
of the boundary conditions. It can be understood by analyz- 
ing the position of rs relative to rsync- We expect that Ci^O 
will have an effect on i™'(,ore oi^ly if '"sync < '"s and the full ef- 
fect will be obtained for rend < '"s- Because rs increases with 
Li,^,coie (grey solid lines in Fig. [1] right panel) and reaches a 
maximum r™' at L"]\^^^, the effect of Ci^O is most prominent 
for Li,^ core close to L™'j,j,j.j.. The right panel of Figure[r|shows 
the effect of CvO on the shock radii. We see that r™' in the 
fiducial calculation (thick black solid line) closely follows the 
results from Paper I (dots) except that the calculations pre- 
sented here have larger r™', because we set dLj^/dr = 0. The 
calculation with CvO (red dashed line) closely follows the 
fiducial results for high M and L^'corc where rsync > '"s"'- 
Here, the Ci^O effect is negligible and i™'core ^^ry close 
to the fiducial value. When r™' ~ rsync, which occurs for 
L^^eore - 1 -5 X lO^^gj-gg g-i ^ 0.32 M© s"'), the CiyO 
effect starts to become important, and both i™*eore and rs 
decrease relative to the fiducial calculation. For L™'core 
2.5 X 10^' ergs s~\ which corresponds to M < 0.035 Mq s~\ 
'"end < '"s"' and CvO affect the structure of the flow signifi- 
cantly and essentially saturate. L™'core with Ci^O is reduced 
to ^-^0.65 of the fiducial value for low M. 

We have seen that the relative positions of rs and 
rsync determine the effect of Ci^O. It is known that 
multi-dimensional effects like convection and SASI con- 
sistently increase shock radii over the corresponding ID 
value and decreas e L™'eore (^ -g- Burrows et al. 1995^, 
Ohnishi et al. 2006; Iwakamii et al.f2008; Miirphy & B urrows 
2008; Marek & Janka 2009; Nordhau££iD itoTol) . In 
Paper I. we parameterized these effects within our steady- 
state models and found that decreased cooling and increased 



heating have very similar consequences except that r™' in- 
creases only in the former case. For this reason, based on 
inspection of simulation results, we argued that the decrease 
of L™'core seen in multi-dimensional simulations is the result 
of less efficient neutrino cooling. Thus, in order to evaluate 
the potentially stronger effect of CuO in multi-dimensional 
simulations, we repeated our calculations, but with the nor- 
malization of the charged-current cooling rate decreased by a 
factor of 2. This decreases i™*core ^y about 40% compared to 
the fiducial case, similar to the difference in i™'core observed 
by iNordhaus et al.l (120101) between ID and 3D simulations. 
The blue dash-dotted line in the left panel of Figure [T] shows 
L™'eore with Ci^O and reduced cooling. As expected, because 
of lower L™'coie and higher rs at fixed L^^^cok, the effect of 
CvO starts to be apparent for M < 0.5 Mq s"' and reaches 
fufl strength forM < 0.04 Mq s"'. 

How does the strength of the Cz/O effect scale with M and 
r„7 In iPaper II we showed that rf^/r^ oc (it^'core)""^^ and 
that the dependencies on other parameters like M and are 
much weaker. We plot values of r™' from IPaper II as dots 
in Figure [T] rig ht pane l, for many different M and r^. We 
also showed in IPaper II that the critical luminosity scales as 
^""core « M° ™M'-^'*r;'-^'. Thus, increasing M or M in- 
creases i™'coie, which in turn decreases r™' and increases rsync 
and rend (eq. IS)) and the effect of CvO will be weaker Sim- 
ilarly, lower r^ yields higher L™'core and lower r"^^ /r^. At 
the same time, r^y^c/ru and r^^^i/r^ will increase (eq. ||6l) and 
the CvO effect will become important at smaller M (smaller 
^™'coie)- Therefore, the collective oscillations will be most 
prominent for the sets of parameters that minimize L"^'-^^^^: 
small M, M, and large r^, . 

In Figure [1] the effect of CvO becomes relevant at M ^ 
0.1 Mq s"'. In fully consistent time-dependent calculations. 
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Figure 2. Relative reduction of i-™'core ^'''^ respect to tlie fiducial ID calculation (/icd) as a function of M including various physical effects. Red dashed fines 
show /led for CvO. Lines are labeled with M and of the PNS. Blue dash-dotted line shows the effect of CvO and reduced cooling rate by a factor of 2 relative 
to a calculation with reduced cooling only. The upper thick grey line shows /red when neutrinos from cooling of the accretion flow are taken into account and the 
lower thick grey fine illustrates /led for cooling rate reduced by a factor 2 ( Paper H) . The grey solid lines wi th points show /ed for multi-dimensional effects when 
the dimension of the simulation is increased from ID to 2D or from ID to 3D (from lNordhaus et aU20l"0l) . 



the values of r^, = 60 km and M = 1 .2Mq chosen for our calcu- 
lation correspond to the very early stages after bounce, when 
M through the shock is still high. For lower r,^ and higher 
M, which correspond to later times after collapse, we plot in 
Figure |2] the ratio of i™'cofe including CvO to i™'core foi" ^ 
reference calculation (/red) - the reduction factor of i™'core 
due to CvO (red lines). We note that in order to evaluate 
only the effect of CvO, we choose the reference calculation to 
have the same M, r^, and microphysics. The blue dash-dotted 
line is obtained for a reference calculation with M = \ .2Mq, 

= 60 km, and cooling reduced by a factor of 2. The two 
thick grey lines have reference calculations without radiation 
transport and with our fiducial neutrino cooling, respectively. 

Figure [U shows that to get a > 10% reduction in i™'core 
for = 60 km and M = 1.2 Mq due to CvO, we require 
M < 0.1 Mq s"'. To get the same reduction in i™'core for 
M = 1.4M0 and = 40km, Figure |2] shows that we require 
M < 0.01 Mq s"'. For = 20km and M = 1.6Mq, there is 
no noticeable reduction in L^'corc because r™' < r^y^^ for the 
whole considered range of M. Looking at solar-metallicity 
supernova progenitor models of IWooslev et al.l (120021^ , M = 
0.1 Mq s"' is reached ^0.65 s after the initiation of the col- 
lapse for an 1 1 .2Mq progenitor, but at ^4 s for a 15 Mq pro- 
genitor. At these times, the accreted masses are M = 1.35 

"^ Ihttp : //www ■ stellarevolution . org/data ■ shtml| 



and 2.2Mq for the 11.2 and 15.0Mq progenitors, respec- 
tively. The lower limit of our calculations M = 0.01 Mq s"' 
is reached only after '-^ 15 s for t he ll.2M pt progenitor, when 
the PNS has almost fully cooled (iPons et al. 1999). From this 
investigation we conclude that the decrease of L™'core '^^^ 1° 
CvO is noticeable only for very low mass progenitors, which 
reach low M at early times, when is still potentially large. 
This would be possible for a stiff equation of state of dense nu- 
clear matter, which would keep ri, higlQ. Multi-dimensional 
effects, which we model with decreased cooling efficiency in 
the flow, would increase rs and thus increase the M where /,ed 
decreases as a result of CvO (blue dash-dotted line in Fig.|2]l. 

Finally, we compare the Ci^O effect to other physical effects 
that have been shown to decrease L^'core- Figure|2]we plot 
with a thick grey line /,ed that was obtained by including the 
luminosity from the cooling of the accretion flow ("a ccretion 
luminosity", dL^/dr=/0) as was calculated inl Paper j We see 
that the decrease of L™'core by 10 to 20% is somewhat smaller 
than the maximum effect from CvO, but is most prominent 
at high M, and th at it is of similar importance at small M and 
small dPaper J) . The lower thick grey line shows /red for the 

^ Note however, that to get an explosion ^"."core ^'^ reached by the 
actual core luminosity L,/,..corc, which depends on the equation of state in 
a more complicated way. Indeed, a softer equation of state generally leads 
to higher Li^,..corc at early times after bounce and may thus be favorable for 
explosion via the neutrino mechanism (e.g. Marek & Janka 2009). 
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cooling rate reduced by a factor of 2, which we use as an ap- 
proximation of multi-dimensional effects. The grey lines with 
dots illustrate the decrease in L™*(,ore due to multi-dimensional 
effects as calculated by Nordhaus et al. (2010). The relative 
drop when going from their ID calculations to 2D is compa- 
rable to the effect of the accretion luminosity, and /red in 3D 
is further decreased. At low M, the effect of CvO becomes 
comparable to that of increasing the dimension of the simu- 
lation from ID to 3D, but only for fairly large and small 
M. 



CM smaller, b ecause physical scale of conversion is moved 
to larger radii (lEsteban-Pretel et al.ll2008t IChakrabortv et all 
lIoTTalbh . 
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4. DISCUSSION & CONCLUSIONS 

We investigate the effect of collective neutrino oscillations 
on the neutrino mechanism of core-collapse supernovae as pa- 
rameterized by the critical neutrino luminosity L™*(,ore- We as- 
sume that the i^g and energies increase to ly^ between a syn- 
ch ronization radiu s rgync and end radius rend, as summarized 
by iDasgupta et al.l (1201 lb . We found that collective oscilla- 
tions affect L™'core if '"sync < rf^, where r"" is the shock radius 
at L™'core' ™d reduce i™'core by a factor of ^ 1 .5 (/red = 0.65, 
Fig. |2|i if Tend < '"s However, this occurs only for low M, 
small M, and large (see Fig.[Tl left panel; Fig.[Tl red lines). 
The reduction of L™'core depends on the assumed energy dif- 
ference between i/^, and Vy,. If the energy e^^ was in- 
creased to 25MeV or 35MeV while keeping e^^ = 13MeV 
and = 15.5 Me V, L^'^ore '^^ ^ would be reduced by 
a factor of 2 and 3.5 (/ed — 0.5 and 0.3), respectively. The 
synchronization radius is rather insensitive to small changes 
in the neutrino energies, because /i(r) cx r~'^ (eq. ID) and thus 
any changes by a constant factor in equation ([T]l have a small 
effect on rsync- 

In order to get a significant decrease of L™'(.ore ^o Ci^O, 
the mass accretion rate needs to be sufficiently low while Vi, 
remains large. This is best achieved in the lowest mass pro- 
genitors, in which M decreases very rapidly due to their steep 
density structure. However, for times < 0.65 s after the col- 
lapse collapse is initiated, the mass accretion rate is likely still 
too high to cause a Ci/O-driven explosion i n an \ \.2M(^ pro- 
genitor (see also IChakraborty et al.ll2011allbt IDasgupta et al.l 
1201 ll) . At late times, M decreases rapidly, but the effect of 
Ci/O on L™'core is likely offset by the simultaneous drop of 
as the PNS cools and the concomitant increase in M. If the 
decrease in is slow enough, perhaps due to a stiff equation 
of state, M might drop enough so that the decrease in i™'(,ore 
due to Ci^O becomes important before black hole formation 
or explosion via the ordinary neutrino mechanism. 

Finally, we note that our implementation of Ci^O is quite 
favorable for having an effect on U^^^^^^^, because we assume 
that the flavor conversion and the increase in the electron neu- 
trino energy occurs for all neutrinos. The flavor conversion 
may not be complete for various reasons. In particular, re- 
finements in the treatment of the neutrino radiation transport 
(multi-group and multi-angle) typically make the effect of 
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